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ABSTRACT 


The results of two different methods for determining the 
rigidity modulus of a soft sediment are compared. In one method 
the resonant characteristics of a torsionally oscillating rod 
which are sensitive to the shear acoustic impedance of the sedi- 
ment in which the rod is imbedded determine the complex rigidity. 
The second method utilizes the observation of the phase velocity 
of an interface wave at the water-sediment boundry. Shear wave 
speeds computed from the experimental data from both methods are 
quite similar in magnitude. For the sediment used here, the 
average value of shear wave speed determined from the interface 
wave sega Pees was 33 m/sec while the shear wave speed determined 
from the measured rigidity was 29 m/sec. The difference lies with- 
in experimental uncertainty. Trends in the mass-physical properties 
of. the sediments are investigated by comparing graphically the de- 
pendence of both the real and imaginary parts of the complex rig:id- 
ity on density, porosity, sound speed, silt and clay percentages, 


Poisson's ratio and density times sound speed squared. 
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I. INTRODUCTION 


Growing concern for the accurate and compr ehensive description 
of the acoustic reflection and attenuation characteristics of the 
ocean bottom is continually stimulated by attempts to fill the 
large gap in predictibility of the actual reflection losses. The 
development of a model for the reflection coefficient which more 
accurately describes the physics of transmission and reflection 
should result in improved capability for predicting transmission 
loss for sound in the ocean in those cases where bottom interaction 
is important, such as passive and active sonars in shallow water 
and active sonars using the bottom bounce mode. 

Bnetorieaie: the ocean bottom has most often been modeled as 
an interface between two fluids of differing densities and sound 
speeds. Usually referred to as the Rayleigh model, it is generally 
accepted because normal surficial sediments have high porosities 
and low rigidities which approximate those of a fluid. The 
Rayleigh model is sometimes found to be deficient, especially at 
lower grazing angles [3,9]. 

Some have accounted for this disparity by presuming that the 
sediment posseses a finite rigidity and hence is capable of sup- 
porting-a shear wave hile Using a viscoelastic model and complex 
Lame' constants, one is capable of describing the formation of 
shear-interface waves and absorption of compressive sound waves 


in the sediment. Bucker, et al. 3] finds better agreement with 








exper imental measurements for the theoretically calculated bottom 
loss using this viscoelastic model and associated sediment para- 
meters obtained by Hamilton Ber koa than for the Rayleigh model 
in a number of cases. 

In order to obtain both complex Lame' constants for the vis- 
coelastic model it is required that, in addition to the measure- 
ment of sound absorption, both shear wave speed and compressional 
wave speed be measured. Hutchins [12] demonstrated the feasibility 
of measuring the complex shear modulus of soft sediments (artifi- 
cial) by noting the effects caused by the rigidity on the resonant 
behavior of a torsionally oscillating rod imbedded in the sediment. 
Utilizing a laboratory prepared kaolin-water mixture, he developed 
the viscoelastometer technique in the frequency range of 38.3 to 
38.9 KHz. 

Cohen Cs], by improving the design of the instrument, increased 
the range of operation to 5.8 to 38 KHz. Cohen used both kaolin- 
water and bent onite-water combinations as artificial sediments. 

Utilizing Hutchins! and Cohen's observations that soft, 
clayey, simulated sediments did, in fact, have a measureable shear 
modulus, Bieda [2] continued this ee earicn by developing a new 
viscoelastometer suitable for use in making in situ measurements 
and the subsequent laboratory analysis of 20 sediment cores ob- 
tained in the Monterey Bay area. 

The results obtained from use of these probes [2,5,12] are 


quite encouraging and are certainly within reasonable bounds for 








- 


this type of sedimentary material (E.L. Hamilton, personal 
communication). However, there has been no independent deter- 
mination of the rigidity for comparison and verification of the 
method. Consequently it was required that experiments be con- 
ducted that would test the validity of the probe method. Suc- 
cessful development of such a probe would make the necessary 
measurements of ocean bottom parameters a more convenient 
operation. 

Another method for determining the elastic constants involves 
measuring directly both the compressional wave and shear wave 
speeds; knowing these, it is a simple matter to determine the real 
elastic constants. The imaginary parts of the elastic moduli must 
be obtained from measurements of the attenuation of the compress- 
ional and shear waves. The former method is used in this experi- 
ment as an independent means for validating the measurement of 
the real components of the rigidity made by the viscoelastometer 
developed by Bieda [2]. The method described is essentially the 
same as that used by Hamilton, et al. [11] in measurement of the 
interface wave speed and compressional wave attenuation in sedi- 


: , i: 
ments in the San Diego area . 


+ Some controversy appears to exist as to just what title should 
be applied to the interface wave that exists at the suface between 
an elastic, solid bottom and the overlying water. Although gener- 
ally referred to as Stoneley waves, it appears that Stoneley 7) 
never included the possibility that the upper medium could be a 
fluid. Other people tend to choose the name Scholte waves since 
Scholte [16] generalized the nomenclature to include fluid-solid 
interface waves. Since the environment in which the author's data 
were taken did not entirely approach an elastic-solid bottom of 
infinite extent below a infinite half-space of homogeneous fluid, 
it was decided to use the term "interface wave'! when describing the 
elastic wave existing at the interface between the bottom and the 
overlying water. 
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II. EQUIPMENT AND PROCEDURES 


A. LOCATION OF THE EXPERIMENT 

Cohen [5] has shown the viscoelastometer to be ineffective in 
sediments of large mean grain Size, possibly because the grain 
dimension becomes comparable to the shear wavelength in the sandy 
sediments and therefore the assumption of homogeneity becomes 
invalid. Thus there was a need for finding a sediment which was 
Silty clay or clayey silt. 

There are vast areas of mud bottoms in the Monterey Submarine 
Canyon area of Monterey Bay, California, that would normally 
suffice for this type of experiment. However, since the services 
of a Pionedsd bile platform were not available the evaluation had 
to be made at a quite shallow depths (1 to 1.5 m) near shore. A 
Suitable mud sediment Site was found in the tidal flats of the 
Kirby Park area of Elkhorn Slough, a few miles north of the city 


of Monterey. 


B. INTERFACE WAVE MEASUREMENT TECHNIQUE 

Direct measurement of the speed of shear waves in very soft 
Marine Sediments is made virtually impossible at acoustic fre- 
quencies because of the rapid attenuation of the shear waves and 
the rere uitey of coupling a suitable transducer to the medium. 
It has, however, been determined by Strict and Ginzbarg [18] and 


by Hamilton, et al. [11] that the ratio of the shear wave speed 
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to the interface wave speed remains relatively constant for a 
Qiven value of the ratio of the two densities of adjoining media 
irrespective of other varying bottom constants. 

While compressional wave speed and wet density of the sediment 
can eaSily be measured in the laboratory and corrected to in situ 
values, it was required that the interface wave speed be measured 
and corrected to shear speed as previously accomplished by 
Hamilton, et al. i aleiele 

An interface wave, generated by the detonation of a blasting 
cap at the water-sediment interface, is examined using three or 
four vertically-polarized, velocity-sensitive geophones mounted at 
known intervals on the interface. An Aluminum rod and flexible 
lines served as a convenient support and the means for spacing of 
the geophones and the electrically-fired blasting caps which were 
mounted under a heavy brass plate. 

During the firing of these blasting caps, the response of the 
geophones is monitered at a remote station. The signal from the 
geophones passes through an impedance-matching transformer, a set 
of attenuator pads and a low-noise amplifier capable of as much as 
50 db gain to a Honeywell Visicorder where the subsequent amplitude 
versus time responses are plotted as shown in Figure 2. A sketch 
of the operational assembly is shown in Figure l. 

The Visicorder traces are studied to determine the points of 
constant phase in each of the geophone traces (point a in Figure 2). 


This is generally best accomplished by attempting to match waveforms 
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which is somewhat aided by irregular spacing of the geophones and 
adjustment of gain so that amplitudes are comparable. Having 
found the points of constant phase, and using convenient time 
reference marks, the interface wave speed is easily determined. 

Having obtained the interface wave speed, the graphs provided 
by Strict and Ginzbarg [18] are entered with the ratio of the 
density of the overlying sea-water to the density of the sediment. 
The graphs are plotted with the ratio of interface wave speed to 
shear wave speed versus the ratio of Shear wave speed to compress- 
ional wave speed, varying the parameter ratio of water density to 
sediment density. Since both the ordinate and abscissa are ratios 
involving the unknown, shear wave speed, an additional assumption 
about the magnitude of shear wave speed allows one to accurately 
determine the necessary value. As the ratio of shear wave speed 
to compressional wave speed approaches 0.5, the individual graphs 
all begin to approach a horizontal asymptote which may conveniently 
be read from the ordinate with a maximum possible error of 
5 percent. 

Having determined the value of the ratio of the interface wave 
speed to shear wave speed, and knowing the interface wave speed, 
it is easy to compute the shear wave speed we = V/A; where A is 
the ratio, ye the shear wave speed and ve the interface wave 


speed). 
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C. THE VISCOELASTOMETER 
1. Equipment 

The apparatus designed by Bieda [2] is a result of a 
series of improvements on the original devices sofatmueene and 
evaluated by Cohen L5] and Hutchins biz). It consists of a 
10%-inch long, *%-inch diameter cylindrical rod which is excited 
in torsional oscillation by a centrally mounted barium titanate 
transducer. The rods, constructed from a constant modulus alloy, 
NI-Span-C, were heattreated and highly polished. The electrical 
leads to the transducer are passed through a supporting yoke which 
Supports a clamp that holds the transducer in position. The clamp 
is secured to the transducer at its center, a node of torsional 
motion for all resonant modes in the rod. A sketch of the 
apparatus is shown in Figure 3. Various compounds were used to 
maintain the water-tight integrity of the electrical connections, 
generally RTV insulating rubber was deposited in the large crevices 
at the junction of the clamp and the transducer while the other 
boundries were coated with varnish or glyptol. 

A Dranetz 100B Impedance-Admittance meter, a suitable 
oscillator and a frequency counter were used to measure the elec- 
trical conductance and the resonant frequencies. 

2. Theory and Method of Measurement 

‘When the viscoelastometer is excited at its mechanical 
resonance, a torsionally oscillating standing wave exists along 
the axis of the rods and transducer. The torsional motion of the 


walls of the rod generates a shear wave which propagates radially 
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outward into the sediment. The radiation impedance of the sur- 
rounding sediment has a modifying effect on the magnitude of the 
electrical xesistance of the probe and its frequency of mechanical 
resonance, compared to those values measured in air (approximately 
no-load conditions). The change in resistance and resonant fre- 
quency are thus direct functions of the rigidity of the sediment [14]. 
Under the assumption that the shear waves that propagate 
outward from the viscoelastometer have wavelengths much smaller than 
the radius of curvature of the rods and that the amplitude of these 
shear waves 1S very rapidly attenuated, the waves propagated in the 


sediment may be considered to be plane waves with the formula: 


T = T, exp [ Yrt/7 (1+i)z J (1) 


where T is the shear stress amplitude, IT, the initial amplitude, 
f the frequency, and nm the shear viscosity. The harmonic time- 
varying function has been suppressed for sake of convenience. 

In a Newtonian fluid, the coefficient of shear viscosity 
is 7 = r/S, where S is the shear strain and S = 0S/dOt. In our 


viscosity model a complex rigidity (G.) has been assumed along 


with a complex viscosity (1.)- Thus, 


Nn. =1, - in, andG_=G, + iG,. (2) 


In a Newtonian fluid, 7, = 0, and 7, = 0 in a perfectly elastic 
solid. Since 7). =T/iWS = - iG / W for simple harmonic motion, 
where Wis the angular velocity, 


G, = &, and G, = wn, . (33) 
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The specific radiation impedance (Z) of the sediment is 
the ratio of shear stress to particle velocity. Following 


Mason L14] and McSkimin C15], for a Newtonian fluid 


Z=R+ ix = rfp (1 + i) (4) 


where R and X are the specific acoustic impedance and reactance 
of the sediment, respectively, and P is the sediment density. For 


the case when viscosity 1S complex, the following relations result: 


_ 2RX _ R°-x* 
"1 ~ @o No ~ “Oo 
(5) 
— R~-x* Be, Cae 
1 p Se oe 


Thus, it becomes apparent that determination of the complex 
rigidity or viscosity requires only a knowledge of real and 
imaginary parts of the specific acoustic impedance for plane 
Shear waves. 

Mason [14] has shown that constants of proportionality can 
be determined which relate the moral ceouct ie resistance’ (R,) to 
the change of electrical resistance (AR); 5 re AR, = K,Ry> and 
the total acoustic reactance (X,) to the change of resonant fre- 
quency (Af), 1.e., Af = KX,» when the loading on the rod is 
changed from air to the medium under study. Here, R = R,/1 and 


xX = X,/1; where 1 is the length of the probe immersed in the 


sediment. Thus, measurement of the acoustic impedance is reduced 





2 : : 
Transducer which drives the rod. 
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to measurement of the change observed in the frequency of 
resonance and in the resonant value of electrical conductance 
when the probe is operated in the material to be analyzed. 
Evaluation of the constants Ky and K, is aeons by 
Measurement of AR, and Af in Newtonian calibration fluids. Since 
uP = O in a Newtonian fluid, R = X; knowing the viscosity and 
density of the test fluid, the procedure for calculating K. and 


ib 


K, is straightforward. 
3. Core Analysis 

In the same relative location of the interface wave experi- 
ment ten ccres were taken for purposes of viscoelastometer analysis, 
density and compressional wave speed measurement. While maintaining 
the core at its original temperature, the complex rigidity is de- 
termined using Bieda's equipment. The measured real part of rigid- 
ity and the observed density are used to determine a shear wave 

L 
speed, i.e., \ em (G,/p)*. 

The core samples were analyzed at two of the resonant 
frequencies, 27 kHz and 38 kHz, due to anomalous behavior at the 
lower frequencies, especially in the stiffer samples. Since the 
constants K, and Ky tend to vary with frequency and viscosity, an 
average value was used at each pay ee from calibrations in 
fluids with viscosities roughly equivalent to those expected in 


the sediments. Table III contains examples of the variance of 


these constants over temperature and viscosity. 
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The computed shear velocities are averaged and are con- 
pared to the value for shear velocity obtained in the interface 


wave experiment in Table II. 


D. COMPRESSIONAL WAVE SPEED MEASUREMENT 

Compressional wave speed is measured in each of the sample 
cores by determination of the time required for a 2 MHz con- 
pressional pulse to pass through a known length of sediment. 
Barium titanate disk transducers are used as a transmitter and 
as a receiver, while the transmitted pulse width and frequency 
are varied in close proximity to 2 MHz in order to obtain resonant 
operation of the transducers. 

The two transducers are mounted facing each other as shown 
in Figure 4. The sample, approximately 3 inches long, is inserted 
in the holder which is in contact with one transducer; the other 
transducer can then be moved to contact the other end of the 
sample. A dial indicator is used to measure the separation while 
a time mark generator and the output of a triggered pulsed oscillator 
are displayed on the face of a dual trace oscilloscope. The tem- 
perature is maintained at a constant value, identical to that 
measured during the interface wave experiment. 

Calibration of the velocimeter is accomplished by measuring 
a sequence of time delays for various transducer spacings ina 
mixture of 22 parts ethanol to 100 parts distilled water. This 
combination was shown to be relatively insensitive to small tem- 
perature changes in the neighborhood of the temperatures of this 


experiment (f2aee 
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The sediment sound speeds and the speed of sound in the water 
samples (computed using Wilson's equation) taken in the vicinity 


of the experiments are shown in Table I. 


E. WET DENSITY AND POROSITY 

Wet density is measured by weighing a known volume of the 
saturated sediment. The weight of this same sample after heating 
in an oven at 105°C for a period of 24 hours is again determined. 
If it is assumed that the weight lost is completely derived from 
evaporation of the entrapped water and that the weight of 1 mé of 


water is lg, the porosity (P) is easily calculated, 


Volume of Voids 
= NORE reE eles. 6 
sf Total Volume cS) 
F, CALCULATED PARAMETERS 
The shear wave speed is calculated by use of the real part of 


the measured rigidity, G and the sediment density, OP, in the 


sad 
formula: 


a a oe (7) 


Bear wave speed is also computed from inteface wave speéd by the 
formula: 

eng V,7A : (8) 
where A 1s the constant of proportionality derived from the graphs 
of Strict and Ginzbarg [18]. Poisson's ratio 0, 1s calculated by 


the relation: 


yds 1 
OS || he (9) 


2 
ONT) -1) 
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where Ms is the compressional wave speed and V. is the shear wave 
speed. The measured and computed physical constants are shown 


im Table [. 


G. SAND-SILT-CLAY RATIOS: MINERALOGY 

Sand-silt-clay ratios were calculated at the Naval Postgraduate 
School by the sieve pipette method. Clay mineralogy analyses, 
performed by Dr. W.R. Bryant, Texas A&M University, are shown in 


the Appendix. 
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III. LIMITATION OF THE METHODS 


A. RESONANT VISCOELASTOMETER TECHNIQUE 

While usage of this probe is not strictly limited to the 
laboratory, a very strong dependence on temperature of the con- 
stants K, and K, (Table III) requires that measurements for com- 
parison be made in a temperature-controlled environment. To 
allow for correct measurements to be made from 2-3/4 inch diameter 
cores, the assumption is made that the shear wave length is much 
smaller than the distance from the probe to the core liner. It 
was noted that at the lower frequency, especially in stiffer 
sediments, erratic readings were generally obtained, possibly from 
the violation of this assumption. This limitation should not pose 
any great problems when the probe is used in situ as there should 
then be no surrounding barrier in coupling to the motion of the 
oscillating rods. . 


It was also noted that the values of Ky and K,, vary with 
changing viscosity as well as frequency and temperature. Con- 
sequently, the final constants used were obtained from an average 
of those values obtained over the lower values of viscosity 
(1-15 poise). This could cause an error in K, of no more than 
30 percent and perhaps 110 percent in the value of K,- 

Another obvious limitation to the laboratory analysis is that 
no matter how carefully a core is taken, there is always some 


structural distortion of the sample as well as some drying before 


analysis can be made. 
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B. INTERFACE WAVE TECHNIQUE 

The most serious breach of assumptions occurs in the measure- 
ment of the interface waves at shallow depths. The determination 
of the interface wave speeds are all predicted on an interface 
between liquid and elastic-solid half spaces ale In the mea- 
surement of the lower frequency (about 10 Hz) interface waves in 
shallow water it would seem that the effect of the overlying water 
1s quite small, since the wavelength for 10 Hz compressional waves 
in water is so large compared to the depth. Thus, the situation 
1S very nearly that normally encountered for Rayleigh waves. This 
effect could cause the ratio of the interface wave speed to the 
Shear wave speed, A, to approach a value of 0.92 (4). This could 
yield a possible error in the measured shear wave speed of as much 
as 10 percent. 

Higher frequency compressional waves, although much faster, 
occasionally tend to mask the arrival of the interface wave due to 
the exponential tail of the explosion. Often this masking caused 
some uncertainty as to the actual arrival time of the interface 
wave.. The standard deviation of the measured values of interface 
wave speed are about 11 percent of the average measured value. 

Since the distance between the geophones enters directly into 
the computation of the interface wave speed, the uncertainty of 
about 2°cm could yield an error of about 5 percent in the result. 

Since only quantitative information describing the relation- 
ship between and interface wave an a solid shear wave is that 


Given by Strict and Ginzbarg [18], the author was forced to make 


Ze 








the approximation that shear wave speed would always be less than 
one-half the compressional wave speed in order to obtain a value 
for the needed ratio, A, as noted in section II-B. This uncer- 
tainty in the constant, A, could yield an error Oe WG greater 


magnitude than 5 percent. 
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IV. RESULTS AND DISCUSSION 


Table II presents a comparison of the two methods of measuring 
rigidity by showing the average values of the shear wave speed as 
calculated from measurements of the real part of the rigidity and 
the shear wave speed computed from the interface wave experiment. 
Extremely good correlation between the two values tendto confirm 
the validity of the torsional viscoelastometer results. [ft 
must be pointed out, however, that the viscoelastometer was 
operated at frequencies of 27 kHz and 38 kHz while the interface 
wave frequencies were about 10 Hz. Although there is a great 
disparity in frequencies, there 1S no reason to expect that 
dispersion of shear speed should occur C11]. The viscoelastometer 
technique is subject to a maximum precision of 30 percent in 
calculating the real part of the rigidity, largely arising from 
an inability to accurately measure all the parameters. While 
these uncertainties of the method point toward the need for 
Stricter maintenance of conditions during viscoelastometer mea- 
surements and during calibration, they are not of sufficient 
magnitude to discredit the method. 

The mechanical disturbance in coring and in transporting the 
Samples .very likely disturbs the physical properties by an inde- 
terminant amount. It is also known from a visual inspection of 
the cores that the bottom in this area was far from being homo- 


geneous, in either the horizontal cr the vertical plane. 
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The measured physical parameters for the ten cores taken tend 
to follow the same general patterns (Figures 5 to 19) as did those 
presented by Bieda [2] with the major exception that the cores 
taken in the slough possess less than 7 percent sand, making 
them much more acceptable for viscoelastometer analysis. 

It is Significant to note that the rigidities obtained ranged 
between 5.1 x 10° BO. 2 571s 10° dynes/em* for the real part and 
ot x 10° to 7.9 x 10° dynes/cem~ for the imaginary part. These 
values, as pointed out by Bieda Lae are significantly lower than 
those found by Hamilton C11], but additional evidence continues 
to point toward the direction of lower rigidities in surficial 
sediments [1,6]. 

The trends noted by Bieda were also observed in this research; 
the real part of rigidity increases with increasing density, 
increasing percent silt and with increasing PV - Rigidity de- 
creases with increasing porosity and percent clay. The same 
trends are noticeable pat less marked with the imaginary part due 


to the greater degree of scatter caused by experimental imprecision. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The relatively close agreement of the two values of shear 
wave speed strongly tends to validate the usefulness of the 
viscoelastometer technique for measuring the real part of the 
complex rigidity in soft, fine grained surficial marine sediments. 

The mass-physical trends noted by Bieda (2) are seemingly 
substantiated by these observations but require a great deal more 
study to completely classify definite behavioral functions. 

In order to further exploit the use of the viscoelastometer, 
it is recommended that further efforts be applied to the design 
of the viscoelastometer in order to minimize the effects on its 
calibration of grain size, stiffness and temperature. 

Further effort is also indicated toward improving the precision 
in measuring the imaginary part of the complex rigidity and toward 
its validation, utilizing a measurement of the sound absorption in 
the sediment as a test. 

Construction of an artificially prepared ocean bottom (kaolin- 
water) is recommended in order to circumvent the damage done to the 
cores when taken from the bottom and to improve homogeneity and 


temperature stability of the medium. 
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TABLE I. Tabulated Physical Properties (cont.) 


G-1 G-2 
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TABLE I. 
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NUMBER 
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Tabulated Physical Properties (cont.) 
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PERCENT 
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TABLE I. Tabulated Physical Properties (cont.) 
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TABLE II. Comparison of Shear Speed Results 


AVERAGE STANDARD DEVIATION ESTIMATE OF 
VALUE OF DATA POINTS PRECISION 

Ratio of V,/V. 0.758 0.062 +5 % 
Shear wave speed 29.0 m/sec 8 m/sec + 30% 
(from measured 
rigidities) 
Interface wave speed 28.8 m/sec 4 m/sec + 10% 
(from interface wave 
experiment ) 
Shear wave speed 33.3 m/sec 4 m/sec + 15% 


(from above) 


Gal 





TABLE III. 


Variation of Constants After Averaging at Low Viscosities 
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Figure 1. Sketch of Interface Wave Measurement Apparatus 
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Figure 2. Visicorder Prints from Interface Wave 
Measurement 
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Figure 3. Sketch of Viscoelastometer and Associated 
Equipment [2] 
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Figure 6. G, as a Function of Porosity 
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Figure 7. G, as a Function of Sound Speed 
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Figure 8. G, as a Function of Poisson's Ratio 
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Figure 9. G, as a Function of Density and Sound Speed 
Squared 
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Figure 10. G) as a Function of Percent Silt 
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Figure 11. G, as a Function of Percent Clay 
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Figure 12. G, as a Function of Density 
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Figure 13. G, as a Function of Porosity 
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Figure 14. G, as a Function of Sound Speed 
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Figure 15. G, as a Function of Poisson's Ratio 
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Figure 17. G, as a Function of Percent Silt 
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Figure 18. G, as a Function of Percent Clay 
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Figure 19. Shear Speed as a Function of Sound Speed 
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APPENDIX: CLAY MINERALOGY 


The results of the 
from X-ray diffraction 
percentage of chlorite 


from areas of outcrops 


quantitative analyses of the clay minerals 
are shown in Table IV. The relatively high 
in these samples is due to river drainage 


of the Franciscan Formation. The graphs of 


G. and G, as a function of the various clay minerals show the data 


ul 2 


from this report as well as those from Bieda (1970)(Figures 20 - 25). 


Little in the way of definite trends may be seen, possibly due to 


the fact that these data are all from nearshore sediment samples. 
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TABLE IV. Clay Mineralogy 


PERCENT PERCENT PERCENT PERCENT 

SAMPLE MONTMORILLONITE ILLITE KAOLINITE CHLORITE 
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Figure 20. G) as a Function of Percent Montmorillonite 
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